Source Tracking for Sco X-1 



K Hayama\ S D Mohanty\ S Desai^, M Rakhmanov\ 
T Summerscales'^, S Yoshida^ 

iThe University of Texas at Brownsville, Brownsville, TX, 78520, US 
^ Center for Gravitational Wave Physics, the Pennsylvania State University, 
University Park, PA 16802, US 
^ . '^Department of Physics, Andrews University, Berrien Springs, MI 49104, US 

^ ^ ' ^Southeastern Louisiana University, Hammond, LA 70402, US 

E-mail: kazuhlr o . hayamaOllgo . org 



o 

O 



l/^ ' Abstract. Sco X-1, the brightest low mass X-ray binary, is likely to be a source 

^S) ' for gravitational wave emission. In one mechanism, emission of a gravitational 

wave arrests the increase in spin frequency due to the accretion torque in a low 
mass X-ray binary. Since the gravitational waveform is unknown, a detection 
^ , method assuming no apriori knowledge of the signal is preferable. In this paper, 

O i we propose to search for a gravitational wave from Sco X-1 using a source tracking 

method based on a coherent network analysis. In the method, we combine data 
KT\' from several interferometric gravitational wave detectors taking into account of the 

direction to Sco X-1, and reconstruct two polarization waveforms at the location 
of Sco X-1 in the sky as Sco X-1 is moving. The source tracking method opens up 
the possibility of searching for a wide variety of signals. We perform Monte Carlo 
simulations and show results for bursts, modeled, short duration periodic sources 
using a simple excess power and a matched filter method on the reconstructed 
^^, signals. 



t^^ ■ PACS numbers: 04.80.Nn,07.05.Kf,95.55.Ym,95.85.Nv,95.85.Sz 
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1. Introduction 



Measured spin frequencies of neutron stars in Low Mass X-ray Binaries (LMXBs) show 

k> , that there is a limiting mechanism on the increase in spin frequency due to accretion 

^ • torque. One possible limiting mechanism is emission of gravitational waves [T]. Then 

5t ! LMXBs such as Sco X-1, which is the strongest X-ray source in the sky, are candidate 

sources of detectable gravitational waves. Also, quasi-periodic oscillations (QPOs) 

have been observed in many LMXBs. Although the mechanism responsible for the 

QPOs is still unknown, there is a model which predicts emission of a gravitational wave 

associated with the QPOs 0. The Rossi X-ray timing explorer (RXTE) satellite has 

millisecond resolution which enables us to probe the physics of transients on rapid 

timescales [J,. It observed QPOs in Sco X-1 at a time LIGO was operational. 

Most current models of these events assume that the gravitational wave from Sco 
X-1 is a continuous periodic signal. Searches for a continuous periodic gravitational 
wave from Sco X-1 were conducted with LIGO S2 data [3], and S4 data [S]. The 
two searches used different approaches: a coherent matched filter method [3] and a 
cross-correlation method [5] . It is also possible that the nature of a gravitational wave 
from Sco X-1 is different from a continuous periodic gravitational wave. 
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In this paper, we adopt a source tracking approach for unniodeled sources based 
on the coherent network analysis RIDGE [^ and propose to monitor Sco X-1 for possible 
gravitational wave emission. In RIDGE, data from several gravitational wave detectors 
are combined coherently, taking into account the antenna patterns, locations of the 
detectors and the sky direction to the source. RIDGE reconstructs h+ and hx time series 
by solving an inverse problem of the response matrix of a detector network. Then we 
do source tracking by reconstructing the h+ and hx time series at the location of 
Sco X-1 in the sky as Sco X-1 is moving. We then apply different algorithms to 
the reconstructed /i+ and hx time series. We perform Monte Carlo simulations and 
show results on the detectability of unmodeled bursts and long, quasi-monochromatic 
gravitational wave signals from Sco X-1, using an excess power and a template based 
search. 

2. Implementation of the source tracking method 

We implement the source tracking of Sco X-1 using a regularized coherent network 
method RIDGE described in [6]. The RIDGE pipeline consists of two main components: 
data conditioning and the generation of detection statistics. The aim of the 
data conditioning is to whiten the data to remove frequency dependence and any 
instrumental artifacts in the data. In RIDGE, the data are whitened by estimating the 
noise floor using a running median [J. By using a running median, the estimation 
of the noise floor can avoid the influence of a large outlier due to a strong sinusoidal 
signal. The whitening filters are implemented using digital finite impulse response 
filters with the transfer function |r(/)| = l/^/S{f), where S{f) is the power spectral 
density (PSD) of detector noise as a function of frequency /. The filter coefficients are 
obtained by using the PSD obtained from a user-specified training data segment. The 
digital filter is then applied to a longer stretch of data. Narrow-band noise artifacts 
known as lines are a typical feature of the output of interferometric gravitational wave 
detectors. Lines originate from effects associated with the functioning of the detectors 
such as mirror/suspension resonant modes and power line interference, calibration 
lines, vacuum pumps etc [5]. These lines dominate over a specific band, which reduces 
the signal-to-noise ratio and the accuracy in recovering gravitational wave signals. In 
RIDGE, the whitening step above is followed by a line estimation and removal method 
by using a technique called median based line tracker (MBLT, for short) described in 
[9]. Essentially, the MBLT method consists of estimating the amplitude and phase 
modulation of a line feature at a given carrier frequency. The advantage of this 
method is that the line removal does not affect the burst signals in any significant 
way. This is an inbuilt feature of MBLT which uses the running median for estimation 
of the line amplitude and phase functions. Transient signals appear as outliers in 
the amplitude/phase time series and are rejected by the running median estimate. 
The resulting conditioned data is passed on to the next step, which consists of the 
generation of a detection statistic. The basic algorithm implemented for coherent 
network analysis in RIDGE is Tikhonov regularized maximum likelihood, described 
in [To] . Recent studies [HI [TOl [12] show the inverse problem of a response matrix of a 
detector network becomes an ill-posed one and the resulting variance of the solution 
become large. This comes from the rank deficiency of the detector response matrix. 
The amount of rank deficiency depends on the sky location, and therefore, plus or 
cross-polarized gravitational wave signals from some directions on the sky become too 
noisy. In RIDGE, we reduce this ill-posed problem by applying the Tikhonov regulator 
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Figure 1. Averaged antenna pattern _Fav to the direction of the Sco X-1. Zero 
time corresponds to 11:00am (Sep. 12, 2007 (UTC)). HI, H2 are the 4 km, 2 km 
LIGO Hanford detector with the same location and orientation, LI is the LIGO 
Livingston detector, and VI is the VIRGO detector. 



which is a function of the sky location. The input to the algorithm is a set of equal 
length, conditioned data segments from the detectors in a given network. The output is 
h+ and hx time series at the sky location of Sco X-1 in the sky, which are reconstructed 
as Sco X-1 is moving. Concatenating the segments of the reconstructed /i-|- and hx time 
series, one can monitor Sco X-1 for bursts or other types of gravitational wave signals. 
In the /i+ and hx time series, signals from the direction of interest are enhanced 
compared with ones from other parts of the sky due to the application of the principle 
of maximum likelihood. 

3. Simulations 



3.1. Information about Sco X-1 

Sco X-1 is the strongest X-ray emitting LMXB on the sky. The right ascension (RA) is 
16h 19m 55.085s, dechnation (DEC) is -15° 38'24.9", and the distance from the earth 
is 2.8 ± 0.3 kpc 13J. Figure [T] shows the averaged detector antenna pattern (i^av) of 
LIGO and VIRGO to the location of Sco X-1 as a function of time since GPS time 
873630000 (11:00am, Sep. 12, 2007 (UTC)). i^av is defined as Fl + F^, where F+ 
and Fx are the antenna pattern function of a plus and cross polarized gravitational 
wave. From this figure, for the LIGO only network, the value of the squared antenna 
pattern exceeds 0.5 for about 33% of the time during a day. If VIRGO is added to 
this network, this network coverage increases to 63%, and the detection efficiency is 
much improved [13] ■ Thus the LIGO- VIRGO network is quite effective to monitoring 
Sco X-1 compared with the LIGO only network. 
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Figure 2. The normalized sensitivity and root-sum-squarc of h of the injected 
sine Gaussian signals [Q = 9) used in the Monte Carlo simulations. The open 
circles are the injected signals of h^ss = [2.1,3.5] X 10~^^Hz~^" with the center 
frequencies 235 Hz. 



3.2. Detection of unmodeled bursts 

We performed Monte Carlo simulations to estimate the detection efficiency for Sco 
X-1. The network consisted of the 4 km and 2 km LIGO Hanford (for short HI, 
H2), LIGO Livingston (LI) and VIRGO (VI). For the detector noise amplitude 
spectral densities, we used the design sensitivity curves for the LIGO and VIRGO 
detectors as given in [15] [16] and we kept the locations and orientations the same 
as the real detectors. For H2, the sensitivity is \/2 times less than HI. Gaussian, 
stationary noise was generated (~ 1000 seconds) by first generating 4 independent 
realizations of white noise and then passing them through finite impulse response 
(FIR) filters having transfer functions that approximately match the design curves. 
To simulate instrumental artifacts, we added sinusoids with large amplitudes at 
(54, 60, 120, 180, 344, 349, 407, 1051) Hz. 

Signals with a fixed amplitude were added to the simulated noise at regular 
intervals. The injected signals corresponded to a single source located at RA = 
16.3 hours and DEC = —15.8 degrees, which was approximately the sky location 
of Sco X-1. We assumed that h+{QJc,t) = Aexp{-{2T:fct)^/2Q^)sm{27:fct) and 
hxiQ,fc,t) = Aexp(— (27r/ci)'^/2Q^)cos(27r/ct), where t is time, Q = 9, was the Q- 
value and fc was the central frequency. The signal strength, A, was specified in terms 

of root-sum-square defined as ft-rss = /-oo*^^ ('*+(^) + ''■x(^)) • In this simulation 

we took hiss — [2.1,3.5] x lO^^^Hz^^". Figure H] shows the averaged sensitivity, 
defined as the amplitude spectral density of simulated HI, H2, LI, VI data divided 

1/2 

by -Fav , and the injected signals. Open circles are the injected signals. The x-axis for 
the circles represents the center frequencies and the y-axis represents the /irss of the 
signals. 
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Figure 3. Block diagram for an excess power method on reconstructed h^ 
and hx time series. The reconstructed h+ and hx time series are contaminated 
with frequency dependent noise. They are passed on to the step of whitening. 
The filter coefficients are estimated using 2 seconds of data which do not contain 
injected signals. The whitened, reconstructed /i_|_ {h\) and hx (f^x) are divided 
into 20 ms segments. Then by calculating power, one obtains detection statistics. 
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Figure 4. Receiver Operating Characteristic. Signal strength in group A is 
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lines correspond to a radial distance statistic, and red lines a excess power statistic. 



The resulting simulated data was passed on to the RIDGE pipeline, and hj^ and 
hx were reconstructed. We ran an excess power algorithm [T^ on the reconstructed 
/i+ and hx- The excess power algorithm is a method to find data segments for 
which the power exceeds a given threshold. Figure [3] shows a block diagram for 
the excess power method. The reconstructed /i+ and hx are contaminated with 
frequency dependent noise originating from the detector. First they were whitened 
using the whitening filter implemented in RIDGE. Filter coefficients were estimated 
using 2 second data intervals which did not contain the injected signals. The whitened 
reconstructed signals Wj^ and ft.^ were divided into segments 20 ms in length. Power, 
which was defined as the summation over the squared amplitude of each segment, was 
calculated. Detection candidates were selected as segments with power beyond a given 
power threshold. We calculated a receiver operating characteristic (ROC) by changing 
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Figure 5. Block diagram for a matched filter method on reconstructed hj^ 
and hx- The reconstructed h^ and hx are passed on a whitening step. The 
filter coefficients are estimated using 2 seconds of data which does not contain 
injected signals. The template signal is also whitened by filters with the same 
filter coefficients. A matched filter method is then applied to the whitened, 
reconstructed signals. 



values of the threshold. Figure [4] shows the ROC curves for the injected signals with 
h-css — [2.1, 3.5] X lO^^^Hz" ' . For comparison, we calculated a detection statistic 
radial distance obtained as follows: In RIDGE we calculated a value of the likelihood of 
the data maximized over all possible /i+ and hx waveforms with durations less than 
or equal to the data segments. The maximum likelihood values were obtained as a 
function of and (f> - this two dimensional output is called a sky-map. Using the 
entire sky-map, we constructed a detection statistic radial distance which scales the 
maximum likelihood by the mean location of the same quantities in the absence of 
a signal. ROC curves were included which were calculated using the radial distance. 
ROC curves for the radial distance statistic were calculated using segments of 0.5 sec 
in length. This result shows even a simple excess power statistic can achieve good 
performance. 



3.3. Detection of modeled bursts 

We considered the detection of well-modeled bursts using a matched filter method. 
We used sine Gaussian signals {Q — 9) with the center frequency of 235 Hz, h^ss 
of 7 X 10~^^Hz~^" as a well-modeled burst. Figure [5] shows a block diagram of 
the matched filter method we used. The reconstructed /i_|_ and h^ were whitened 
with the whitening filter. The filter parameters were estimated using the preceding 
2 seconds of data which did not contain signals. The template was also whitened by 
the same whitening filter and filter parameters. The whitened data was passed on to 

the matched filtering defined as T = x^ s"^ / av s 



st^st, 



where x ^ hi {* = [- 



]) and 

s^ were column vectors of the whitened reconstructed polarized waveforms and the 
whitened template, a was the standard deviation of the whitened data. T was the 
signal-to-noise ratio (SNR). 

Figure ?? is the result. The top plots show the reconstructed /i+ (left) and h^ 
(right). The bottom plots are the SNR. These plots show the signals are detected 
with SNR ~ 15, 15 for ft,+ and h^ . Considering SNRs of the injected signal in HI, 
H2, LI, VI are 15.2, 10.7, 18.9, 18.5 respectively, this result is reasonable. 



5.^. Detection of a monochromatic signal 

Finally we demonstrated the detection of a short duration monochromatic signal. 
The frequency of the injected signal was 600 Hz and the /ij-ss of the injected signal 
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Figure 6. Top plots show the reconstructed h^ (left) and hx (right). The 
corresponding injected signal is a sine Gaussian signal {Q = 9) with the center 
frequency of 235 Hz, h^ss = 3.5 X 10~^^Hz~^". Bottom plots are output of the 
matched filter method. 
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Figure 7. Top plots show the reconstructed h+ (left) and hx (right), and 
red lines show injected h^ and hx, oi which frequency is 600 Hz and hrss is 
2.1 X 10~'^^Hz~^''^. Bottom plots show the amplitude spectral density (ASD) of 
the reconstructed /i_|_ (left) and hx (right). 
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was 2.1 X lO^^^Hz^^". The top plots in figure [7] are of the whitened reconstructed 
/i+ (left) and hx (right). The blue lines are reconstructed /i+ (left) and h^ (right) 
and the red plots are the original injected signals. The lower plots are of amplitude 
spectral densities. The reconstructed h+ and hx oi the signal can be seen clearly at 
600 Hz in the plots. 

4. Conclusion 

We proposed the source tracking for Sco X-1 using the regularized coherent network 
analysis method RIDGE. To demonstrate its performance, We considered the detection 
of unmodeled bursts, modeled bursts, and continuous sources. We applied an excess 
power method to the reconstructed h^ and hx for the detection of the unmodeled 
bursts. We also applied a matched filter method to the reconstructed waveforms for 
modeled signals. Finally, we reconstructed the short duration periodic signal, and 
applied Fourier transform to the reconstructed waveforms. These results show the 
source tracking approach works for a variety of signals. 
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